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ABSTRACT 


A  rcviaw  of  the  technical  research  on  composite  Materials 
used  in  the  manufacture  of  solid  propellant  missiles  <3i«e~ been— 
N conducted.  Current  design  methods  were  compared  to  the  methods 
given  in  the  literature.  Host  of  the  designers  of  composite 
structures  use  empirical  relationships  based  on  data  compiled 
from  previous  designs,  to  predict  the  structure's  characteris¬ 
tics.  No  universal  approach  to  characterization  of  material 
properties  of  composites  has  been  found.  because  no  one  has 
been  able  to  adequately  characterize  the  material  prjperties  of 
the  constituent  materials  and  their  interrelationships.  Recent 
analyses  have  shown  that  only  finite  element  analysis  CFEA)  can 
adequately  represent  the  material  loading  within  the  composite 
matrix.  Accurate  material  properties,  obtained  only  by  detailed 
analysis  of  experimental  work.  are  required  to  be  used  in 
the  calculations.  Tha  properties  are  deduced  from  the  experimen¬ 
tal  results  by  computer  analysis  by  finite  elements.  Accurate 
aralyaia  of  the  total  structure  of  a  composite  material  chamber 
will  also  require  the  use  of  3-dimensional  FEA.  The  use  of  3-D 
FE  has  been  limited  because  of  computer  tine  and  coats.  Howeysr. 
new  computers  will  now  handle  reasonably  sized  3-diaensiona.l 
problems  within  acceptable  limits  of  time  and  coat. 
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NOMENCLATURE 


E  Young's  modulus 

G  Shear  Modulus 

K  stress  concentration  factor 

M  moment 

P  pressure 

x,y,z  principle  directions,  cartesion  coordinates 

r,0  principle  directions,  polar  coordinates 

V  volume 

W  weight 

a  stress,  normal 

t  stress,  shear 

e  strain 

6  angular  displacement,  in  degrees 

v  Poisson's  ratio 
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INTRODUCTION 


Composite  materials  have  been  used  for  solid  rocket  motor 
applications  for  almost  three  decades.  The  development  of  these 
materials  progressed  at  a  rate  significantly  ahead  of  the  basic 
theoretical  analysis,  so  a  trial  and  error  development  program 
was  the  foundation  for  years  of  manufacturing.  The  design 
parameters  of  isotropic  materials  are  highly  developed,  and  data 
on  the  material  properties,  including  most  failure  properties, 
are  readily  available.  However,  no  similar  information  is 
available  for  composite  materials.  Designers  must  rely  on  the 
meager  collection  of  data  that  has  been  compiled  from  experi¬ 
mental  programs  with  similar  designs  and  composite  structures. 

Significant  work  has  been  performed  in  the  development 
programs  of  rocket  chambers,  to  test  and  evaluate  sub-scale 
pressure  vessels  of  similar  design  ahd  composite  structure. 
These  sub-scale  results  are  then  scaled-up  to  the  design  require¬ 
ments.  Unfortunately,  the  scaling  factors  are  not  well  under¬ 
stood  and  failures  often  occur;  this  is  primarily  the  result  of 
the  lack  of  fundamental  knowledge  of  the  composite  structures 
material  interactions. 

Study  has  been  devoted  to  this  fundamental  deficiency. 
Researchers  have  attempted  to  apply  basic  isotropic  material 
relationships  to  specific  composite  structures  and  have  had 
limited  success  in  the  determination  of  the  characteristics. 
Isotropic  materials  such  as  metals,  crystals  and  some  plastics, 
are  not  purely  isotropic;  their  failure  properties  account  for 
the  irregularities  and  characterize  the  failure  mode<e>  that 
result  according  to  the  size,  location  and  orientation  of  these 
irregularities.  Since  composite  materials  are,  by  their  very 
nature,  anisotropic,  the  attempt  to  characterize  their  material 
properties  with  relationships  that  relate  the  defects  in  their 
structure-  to  their  strength  and  failure  modes  is  misguided. 

In  this  paper  I  will  survey  the  theoretical  development  of 
composite  structures  and  contrast  the  theory  with  the  current 
manufacturing  design  practices.  I  will  discuss  the  proofing 
methods  end  those  methods  that  are  used  to  analyze  damage  to 
pressure  vessels.  Conclusions  and  recommendations  will  be 
offered  to  improve  the  design  and  analysis  methods  used  in 
current  chamber  design  and  structural  analyses. 


THEORETICAL  DEVELOPHENT/LITERATURE  REVIEW 


Laminate  Theory  <LT) 

The  general  laminate  theory  is  a  Mathematical  means  of 
calculating  the  material  properties  of  a  composite  from  the 
angles  of  fiber  orientation  and  the  properties  of  the  layers  of 
the  fibers.  The  matrix  properties  are  included  in  this  calcula¬ 
tion  and  the  matrix  is  assumed  to  be  homogeneous  and  perfectly 
bonded  to  the  fibers.  The  elastic  properties  of  the  layers  are 
determined  from  the  Young's  moduli  En,  Ei2»  the  shear  modulus 
Gi2»  and  the  Poisson's  ratio  v\2  using  Hooke's-law  relationship: 

o  sEe,  Eqn. 1 

where  strain  is  proportional  to  stress. 

The  laminate  properties  are  calculated  by  using  the  angle  of 
the  laminate  direction  and  the  angle  of  the  material,  which  are 
then  added  vectorially  (see  Appendix  1). 

Theory  implies  that  the  stacking  sequence  and  boundaries 
will  have  no  effect  bn  laminate  strength.  Pagano  and  Pipes  £13 
have  shown  that  this  is  not  true.  •  Their  work  illustrated 
significant  stacking  sequence  effects,  as  well  as  significant 
boundary  effects.  Tolbert  £23  also  reported  that  the  laminate 
theory  failed  in  the  analysis  of  thick-walled  cylinders,  but  was 
able  to  represent  thin-walled  vessels  without  serious  error.  The 
primary  fault  of  this  theory  is  in  its  assumption  of  zero 
boundary  stress  fields.  It  seems  to  accurately  account  for  the 
stress  fields  remote  from  the  boundary  layers. 

Many  researchers  have  attempted  bo  confirm  this  theoretical 
development;  however,,  there  has  been  little  success  in  correlat¬ 
ing  the  experimental  results  to  the  theory.  The  assumptions  made 
to  simplify  the  calculations  do  not  allow  for  the  physical 
defects  in  a  manufactured  laminate  structure  of  a  thickness 
greater  than  one-tenth  of  an  inch  £33 .  , 

For  rough  calculations,  the  general  laminate  theory  can  be 
used  to-  obtain  the  approximate  laminate  strength;  the  result  can 
be  refined  by  experimental  trial  and  error  procedures.  Chamis 
£43  published  a  procedure,  based  on  laminate  theory,  that  can  be 
adapted  for  use  on  small  computers  or  hand-held  calculators. 

The  relationships  that  apply  to  material  properties  are  not 
affected  by  thickness  and  are  discussed  in  the  theory  of  applied 
mechanics.  Fracture  mechanics  change  with  specimen  thickness. 
Fracture  mechanics  theory  is  highly  developed  for  isotropic 
materials,  but  not  for  anisotropic  laminate  structures.  Further, 


the  relationships  between  the  fracture  characteristics  of  thick 
and  thin  fiber-reinforced  laminates  have  not  been  established. 


Linear  Elastic  Fracture  Mechanics  (LEFM) 

Several  investigators  have  modified  LEFM  for  a  homogeneous 
anisotropic  material  (see  Appendix  2) .  They  have  developed 
expressions  for  the  crack  tip  stress  field,  using  Hooke's  law, 
for  a  homogeneous  linear  anisotropic  material  in  the  case  of 
plane  strain  and  pure  shear.  The  results  were  analogous  to  the 
isotropic  material  case.  The  results  for  non-homogeneous 
anisotropic  laminates  were  different,  but  were  characterized  by 
the  "inherant  flaw  model”.  In  this  model,  a  "high  intensity 
energy  region"  adjacent  to  the  holes  and  at  the  tip  of  the 
notches  was  postulated.  This  region  was  assigned  a  dimension 
which  was  then  used  in  the  plastic  zone  correction  factor  for 
isotropic  materials  to  modify  the  stress  intensity  factors  and 
the  stress  distribution.  The  investigators  found  that  single 
cracks  do  not  nucleate  and  grow,  as  in  metals,  under  fatigue 
loading.  This  loads  to  the  question  of  how  applicable  LEFM  is 
for  composite  materials. 

While  trying  to  explain  this  phenomenon  and  the. characteris¬ 
tic  problem,  the  investigators  postulated  a  non-fracture  mech¬ 
anics  stress  fracture  criteria.  A  “point  stress  model"  postula¬ 
ted  that  failure  would  occur  when  the  average  stress  at  some 
distance  from  the  discontinuity  reached  the  critical  stress. 
This  model  worked  well  with  some  laminates,  and  not  so  well  with 
others. 

Using  the  LEFM  and  the  R-curve  concept,  which  relates  the 
rising  fracture  toughness  of  the  material  to  the  increasing  crack 
extension,  the  size  of  the  plastic  zone  for  conventional  mat¬ 
erials  of  thin  sheets  may  be  considered.  For  fiber-reinforced 
thin  sheets,  the  point  and  average  stress  models  explain  the 
increase  in  toughness  as  the  notch  length  increases;  as  a  result 
of  the  stress  distribution  near  the  notch  tip.  These  approaches 
consider  the  deviation  of  the  load/crack  opening  displacement 
curve  from  linearity  because  of  the  damage  zone  formation 
and  self-similar  crack  extension.  This  does  not  apply  in  thick 
laminates  or  in  laminates,  with  delaainations  CS3 . 

Shih  and  Logsdon  C6J  found  during  their  study  of  thick- 
section  graphite/epoxy  composites  that  the  LEFM  did  not  apply 
directly  to  the  composites  which  had  cracks  perpendicular 
to  the  fiber  orientation.  They  also  found  that  in  three-point 
bending  the  failure  mode  was  interply  failure.  This  mode  of 
failure  faults  the  basic  assumptions  in  LEFM. 

Self-simil-  -*  crack-propagation ,  another  assumption  of  the 


homogeneous  anisotropic  mode?',  is  not  assured  because  of  the 
anisotropy  and  heterogev- ai c.y  of  the  laminate  microstructure. 
Only  for  fractures  that  : re  confined  to  interlaminar  planes, 
where  the  matrix  phase  dominates  the  fracture  response  and 
where  self-similar  propagation  is  assured,  can  LEFM  be  used 
effectively  C7] .  In  these  cases,  tests  can  determine  the  stress 
intensity  factors,  Kjj,  and  Km  for  the  opening,  shedring 
and  tearing  modes,  respectively,  as  well  as  the  mixed  mode  of 
fracture. 

Other  researchers  found  that  the  values  of  toughness  reached 
asymptotic  values  and  nonhomogenei ties  because  of  fiber  size 
and  ply  thickness  were  not  significant.  Simple  macromechanical 
models  based  on  LEFM,  such  as  the  boundary  integral  equation 
(BIE)  method,  could  be  employed  to  predict  fractures. 

Further  models,  such  as  the  "shear-lag  model",  were  develop¬ 
ed  to  predict  the  behavior  of  unidirectional  fiber-reinforced 
composites  that  were  susceptible  to  matrix  splitting  ahead  of 
the  crack  tip  vs.  fiber.  This  two-dimensional  model  assumes  that 
the  fibers  support  all  the  axial  force  (because  of  their  higher 
elastic  modulus)  and  the  matrix  supports  the  shear  and  trans¬ 
verse  normal  forces.  This  method  can  be  modified  for  angle  ply 
laminates  also.  The  micromechanical  models,  such  as  the  BIE  and 
the  "shear-lag”  model,  were  generally  more  accurate  in  the 
fracture  predictions  than  in  the  macromechanical  models. 

The  macromechanical  models  which  treat  the  composite  as  a 
homogeneous  anisotropic  medium  do  not  fully  explain  the  physical 
failure  mechanism,  even  when  they  agree  with  the  experimental 
results.  The  micromechanical  models  have  the  potential  to 
account  for  the  effects  of  material  nonhomogenei ty.  These 
nonhomogeneities  can  be  addressed  by  careful  analysis  with  finite 
element  methods  C8) . 


The  finite  element  method  of  failure  characterization 
assumes  that  when  an  individual  finite  element  exhausts  all  of 
its  capacity  for  strain  energy,  it  fails,  forming  a  crack.  The 
element  is  then  removed  from  the  matrix;  the  reaction  loads  it 
had  are  transferred  to  the  adjacent  elements;  its  stiffness 
property  is  reduced  to  zero.  Zf  the  reaction  leads  exceed  the 
capacity  of  adjacent  elementa,  they  too  fall.  Should  the  applied 
stress  on  the  composite  remain  constant  during  this  process, 
catastrophic  failure  occurs.  This  model  can  be  used  to  study  thk> 
basic  failure  processes  that  take  place  in  matrix  material.  The 
model  can  study  post*debond  sliding  between  fiber  and  matrix, 
fiber  fracture,  and  fiber  pullout,  to  mention  a  few  of  the  many 
possibilities. 


Y'^J. 


The  full  characterization  of  the  failure  of  an  inhomogen- 
eoua  anisotropic  laminate  material  is  exceptionally  complex.  The 
planar  characteristics  have  been  formulated  by  Wang  C9J  whoso 
Model  considers  a  partially  closed  delaminntion  having  friction¬ 
less  crack-surface  contact.  The  result  is  the  classical  inverse 
square-root  singularity  with  an  algebraic  multiplicity  of  two. 
With  frictional  contact,  the  results  depend  on  the  frictional 
coefficients  and  the  fiber  orientations  of  the  adjacent  dis¬ 
similar  plies.  If  the  delamination  has  n  very  small  area  of 
crack  closure,  the  model  is  simplified  by  taking  the  limit  of  the 
partially  closed  delamination.  Wang  used  Leknitakii's  complex- 
variable  stress  potentials  in  conjunction  with  an  eigenfunction 
expansion  method,  which  leads  to  a  standard  eigenvalue  problem. 
Frow  this  eigenvalue  problem,  stress  singularities  and  the 
general  solution  structures  of  deformation  and  stress  fields 
associated  with  composite  delamination  are  determined.  The 
complexities  of  the  composite  delamination  phenomenon  include  the 
inherent  crack-tip  singularity,  the  effect  of  anisotropy  of  each 
constituent  fiber  lamina,  and  the  abrupt  change  of  stiffness  or 
ply  orientation  through  the  laminate  thickness  direction.  In 
addition,  the  3-dimensional  state  of  stress  and  deformation  that 
exists  at  the  delamination  gives  rise  to  all  modes  of  fracture, 
vhich  renders  the  problem  mathematically  intractable. 

For  pressure  vessels  and  large  planar  structures,  it  is  not 
only  convenient,  but  accurate  to  model  these  structures  as 
2-dimensional.  This  simplification  allows  accurate  mathematical 
formulation  and  solution  with  finite  element  analysis  and  methods 
as  advanced  recently  by  Wang  MOJ  and  many  other  authors, 

, These  analyses  require  the  material  properties  to  be 
accurate  for  the  fracture  conditions.  This  requirement  will 
force  the  testing  methods  to  accurately  characterize  the  cbndi- 
tions  in  the  composite,  including  the  fiber-resin  interaction  and 
the  resin  properties  in  the  laminate  till  .  The  resin  properties 
in  the  laminate  may  differ  significantly  from  test  coupons 
because  of  processing  variations  and  conditions  that  cannot  be 
duplicated  in  a  test  coupon.  There  may  be  a  aignif icantly  higher 
concentration  of  trapped  air  bubbles  and  alrborna  particlss.  The 
environmental  conditions  and  winding  variables,  such  as  winding 
speed  and  tension,  will  affect  the  resin  properties  and,  hence, 
the  laminate  properties. 


3-dlnecalonal  Finite  Element 


In  the  structures  where  the  atreee  state  i.i  not  planar,  and 
e  true  3-dimensional  stress  state  exists,  the  only  accurate 
analysis  of  the  fracture  behavior  is  by  3-dimensional  finite 
element.  Although  this  method  hea  been  formulated  for  years,  the 
mathematical  difficulties,  because  of  the  sheer  number  of 

3- 


calculations  required,  have  precluded  its  effective  use  until 
recently.  The  structures  of  principal  concern  include  aircraft 
structural  members  and  pressure  vessels  that  are  net  perfect 
spheres.  Cylindrical  pr^uaure  vessels  heve  stress  states  that 
approach  planar  in  whe  cylinder  section,  but  have  complicated 
3-dimensional  stress  states  in  the  dome  end  sections. 

Wang  £123  has  solved  the  problem  of  composite  delamination 
and  failure  using  the  theory  of  anisotropic  laminate  elasticity 
and  interlaminar  fracture  rechanics.  These  solutions  require  the 
use  of  3-dimensional  finite  element  analysis  to  accommodate  the 
3-dimensional  stress  state  and  interactions  of  the  lamintte 
materials.  This  problem  is  very  complex  because  of  the  geometric 
and  material  discontinuities  and  the  inherently  coupled  inodes  of 
fracture.  Wang  presents  the  effects  of  the  geometry,  the 
lamination  and  the  crack  variables.  The  basic  governing  equa¬ 
tions  are  lengthy  and  complex  and  are  well  presented  in  £123 ,  so 
I  omit  them  here.  His  results  do  show  that  the  stress  intensity 
factor  for  the  tearing  mode  is  one  or  two  orders  of  magnitude 
higher  than  the  other  modes.  .  Wang  also  confirmed  that  the  fiber 
orientations,  ply  thickness,  stacking  sequences  and  other 
lamination  variables  all  effect  the  delamination  behavior 
significantly. 

A  more  efficient,  cost-effective  finite  element  method 
was  developed  based  on  the  concept  of  the  mode  I.II  magnitude  by 
Nishioka  and  Altluri  £133 .  Thia  was  a  one-step,  smaller  mesh 
refinement  of  their  hybrid  3-dimensional  method  that  used  a 
two-step  solution  with  a  very  fine  mesh. 


The  sequence  of  the  stacking  of  the  fiber  layers  has  a 

significant  effect  on  the  strength  of  the  composite  £1,103. 

Herakovich  discussed  the  theoretical  implications  of  the  stacking' 
sequence  through  the  theories  analysed  by  finite  element  methods 
£143..  The  use  of  adjacent  t  9  layers  leads  to  high  interlaminar 
shear  stresses  and  that  interlaminar  normal  stresses  are  very 
dependent  on  the  stacking  sequence.  Therefore,  adjacent  ±6 
layers  should  generally  be  avoided.  Interspersed  tQ  layers  are 
preferred  to  reduce  the  extremes  of  interlaminar  snear  and  normal 
streesee.  Thi»  is  I.-scau no  of  a  "mismatch"  of  the  engineering 
properties  batveen  adjacent  .layers.  The  interlaminar  shear 
stresses  are  a  function  of  the  coefficients  of  mutual  influence 
that  can  be  ten  times  greater  then  the  mismatch  of  Poisson's 
ratio.  This  mismatch  reaches  a  peak  at  10°  -  15?  fiber  range  for 
±0  laminae  combinationa.  It  can  be  reducad  by  a  factor  of  two 
whan  the  layers  are  Interspersed  between  8° and  layers. 

Layer  thickness,  both  total  and  fasin/fiber  ply,  affects  the 


laminate  structural  characteristics.  Apart  from  the  obvious 
linear  effect  of  multiple  plies  of  specific  strength.  the 
composite's  fracture  end  failure  characteristics  are  altered  with 
variations  in  thickness-  The  stacking  sequence  of  the  plies  also 
effects  the  overall  thickness;  a  clustered  sequence  of  plies  can 
be  twice  as  thick  as  an  alternated  sequence  as  uhown  in  Figure 
1.  Although  the  in-plane  elastic  properties  of  such  laminates 
are  independent  of  stacking  sequence,  the  strength  and  toughness 
can  vary  significantly  CIS]. 


i  :  . 

(  ■  i 

i  .  I 
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Stacking  Sequence  and  Fiber  Orientation 
Figure  1. 

Resin  systems  are  formulated  tc  provide  specific  matrix 
properties.  processing  characters  .ics  and  fiber 'resin  bonds. 
These  resins  must  have  chemical  characteristics  that  are  compati¬ 
ble  with  the  flber(a)  being  used  in  the  composite.  Damage 
tolerance  of  a  composite  structure  is  a  function  of  the  resin 
volume  content  and  the  plastic  flow  characteristics  of  the 
resin.  This  resin  also  effects  the  compressive  strength  of  the 
composite.  The  failure  mode  frequently  observed  is  a  delamina¬ 
tion  fracture  in  the  through  thickness  direction.  Therefore, 
improving  the  compressive  strength  of  the  resin/fiber  system 
should  improve  the  damage  tolerance  of  the  composite.  Resin 
systems  with  elastomeric  additives,  thermoplastic  additives  and 
viryl  modifiers  have  shown  improved  tolerance  to  damage  [16] . 
The  interlaminar  shear  strength  . of  the  composite  is  dependent 
on  tha  fiber-resin  bond;  fiber  surface  treatment  may  improve  the 
composite's  characteristics  C171  , 


CURRENT  APPLICATIONS 

The  designers  of  the  chambers  for  modern  solid  rocket  motors 
use  many  methods.  The  designer  is  given  the  parameters  and 
allowables  for  a  feasible  design  based  on  previous  data  and 
expected  gains  through  new  materials  and/or  approaches.  The 
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structural  members  of  the  design  are  chosen  to  balance  reli¬ 
ability  and  performance.  The  reliability  is  calculated  by 
the  expected  number  of  failures  under  normal  operating  condi¬ 
tions.  The  performance  of  missile  chambers  is  calculated  from 
the  product  of  the  pressure  inside  the  chamber  and  the  volume  of 
the  chamber  divided  by  the  weight  of  the  chamber  and  is  repre¬ 
sented  by: 

Performance  *  PrV/W  Eqn.2 

The  theories  that  are  applied  to  the  requirements  are  in  two 
categories:  thin-wall  and  thick-wall.  The  typical  missile 
chamber  is  in  the  thin-wall  category,  whereas  the  segments  for 
the  filament-wound  chambers  used  in  the  solid  rocket  motors  for 
the  Space  Shuttle  Program  are  in  the  thick-wall  category.  The 
exact  transition  point  between  the  two  criteria  has  not  been 
fully  characterized. 

With  these  limitations,  the  desi  ler  calculates  the 
loads  the  chamber  will  experience  and,  after  flying  appropriate 
margins  of  safety  for  the  reliability  cri  e^ ion,  designs  the 
required  composite  structure.  He  typically  predicts  the  expected 
fiber  stresses  by  extrapolating  those  that  were  experienced  in 
previous  designs  similar  to  the  one  in  development. 

From  the  data  that  he  has  access  to  on  other  similar 
designs,  the  designer  derives  material  properties.  The  structural 
tasting  of  design  prototypes  during  development  is  crucial.  An 
initial  design  is  then  formulated  by  combining  netting  theory,  a 
winding  model,  and  dome-contour  empirical  equations.  These 
contour  equations  are  varied  through  several  dome  designs  to 
obtain  c  design  that  meets  the  design  criteria.  Two-dimensional 
finite  element  analysis  is  .  applied  to  locate  the  hoop. rsversal 
points,  to  size  port  reinforcements,  and  to  determine  the 
hardware  geometries.  These  finite  element  codes  may  be  separated 
for  the  forward  and  after  dome  regions,  if  they  are  different  by 
virtue  of  the  symmetry  of  the  chamber.  This  separation  may 
significantly  reduce  the  computer  time  required  to  process  the 
model.  The  separation  is  made  at  a  point  on  the  chamber  where 
*he  loading  is  purely,  or  nearly  so,  membrane  load.  .  The  2-D 
finite  element  analysis  provides  the  designer  With  the  displace¬ 
ments,  strains,  and  stresses  for  the  loading  of  the  fiber  matrix 
analysis. 

Safety  factors  and  margins  of  safety  are  applied  in  accord¬ 
ance  with  the  HIL  Handbook  V  Interaction  Equation  based  on  the 
worst-case  <\et  of  peak  loads.  The  skirts  are  designed  for  the 
chamber  after  t initial  design  of  the  pressure  vessel.  The 
skirts  provide  for  the  structural  transmission  of  the  boost 
energy,  as  well  as  providing  the  handling  and  attaching  fixture 
for  the  pressure  vessel.  The  skirts  are  designed  to  resist 
buckling  in  all  modjes  including:  bending,  differential  pressure. 


axial  loading  and  shear. 

Tests  are  required  .  to  adjust  er.pirical  and  sexi-empirical 
formulas  and  relationships  which  will  account  for  designs  &  test 
results  that  differ  from  the  previously  used  values. 

These  adjustments  usually  account  for  the  differences  in 
material  testing  properties  from  standard  tensile  and  coupon 
tests  and  allow  for  the  composite  material  interaction.  Although 
these  adjustments  would  allow  for  the  back-calculation  of 
material  properties  as  they  apply  in  the  composite  structure;  the 
relationships  are  .  usually  adjusted  instea ?  The  adjusted 
relationships  become  very  effective  for  the  design  in  considera¬ 
tion  and  can  result  in  actual  performance  at  failure  being  within 
a  few  percent  of  predictions. 

This  empiricism  is  an  acceptable  and  prudent  method  of 
design  for  those  applications  that  closely  resemble  in  structure, 
or  composition,  a  previous  design  with  a  given  data  base.  For 
radical  designs  or  composite  structures,  the  empiricism  would  be 
as  costly  as  developing  the  theory  to  the  application. 


Three-dimensional  analysis  is  used  when  other  methods  fail 
to  account  for  an  anomalous  condition  or  a  failure  that  pannot  be 
explained  with  existing  analysis.  The  software  and  engineering 
is  available  to  make  this  a  routine  evolution  and  the  basis  for 
most  design  evaluations. 


Current  testing  procedures  include  the  instrumentation  of 
chambers,  real-time  and  standard  radiography,  structural  loading, 
and  the  hydro-proofing  and  hydro-bursting  of  chambers.  The  proof 
testing  of  composite  materials  was  studied  by  Hahn  and  Hwang 
C181 .  The  lifetir*  and  strength  of  most  materials  usually 
show  more  scatter  than  the  other  material  properties,  since  they 

[in  the  homogeneous  structure.  In 


are  related  to  local  defects 


the  composite  structures,  lifetime  is  related  to  static 


strength, 
been  found 

that  can  be  endured  hao  not  basin  exactly  characterized,  although 


The  effect  of  proof 
to  be  negligible. 


no  apparent  reduction  . in  st 
tested  through  the  design  loads 


rength  ha«A  been  noted  in  vessels 


X 


Materials  for  the  composi] 
review.  Current  experimental 
increase  the  performance  of 
order  of  1S-25X.  Laboratory 
fiber /resin  systems  has  had 
results  represent  tests  with 
been  standardized  in  dimensions 
results  would  then  be  scaled  to 
methods. 


testlpg  on  subsequent  strength  has 
The  actual  amount  of  proof  stress 


|te  structures  are  constantly  in 
rk  in  this  field  is  attempting  to 
mposlte  pressure  chambers  on  the 
Resting  of  nee  fibers,  resins  and 
promising  results  (19] .  These 
|small  pressure  bottles  which  have 
for  compilation  of  data.  The 
a  future  design  using  existing 
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CONCLUSION?;  AND  RECOMMENDATIONS 


The  current  level  of  theoretical  development  of  composite 
chambers  is  well  beyond  the  practical  application  level  used  in 
the  design  and  manufacture  of  these. chambers.  The  most  promising 
area  of  fundamental  growth  that  can  be  applied  to  the  designing 
and  manufacturing  areas  is  that  of  accurate  3-dimensional  finite 
element  analysis  <3-DFE).  The  limiting  factors  precluding  this 
transition  include: 

1)  The  reluctance  of  industry  to  adapt  hardware  and 
software  to  perform  these  analyses. 

2)  The  lack  of  expertise  at  the  industry  level  to  ef¬ 
fectively  employ  the  procedures  of  this  method,  including  the 
selection  of  the  regions  where  3-DFE  analysis  is  required. 

3)  The  lack  of  testing  methods  and  systems  to  determine  '.he 
material  properties  within  the  composite  structure  by  deduction 
through  3-DFE  methods;  e.g.  the  deduction  of  the  material 
properties  from  test  results  instead  of  vice  versa. 

The  longitudinal  material  property  values  that  are  compiled 
in  the  literature  are  not  very  useful  in  the  designs;  actual 
designs  include  many  property-reducing  factors  such  as  stress 
concentrations,  fiber  wash,  fiber1  slipping,  resin-rich  and 
resin-lean  regions,  delaminations  and  void*,.  These  factors, 
although  not  part  of  the  design,  exist  in  all  manufactured 
chambers  to  some  degree.  The  analyses  of-  the  current  designs  are 
based  on  the  data  from  previous  tests  and  .  have  empirically 
included  these  effects  on  a  macroscopic  scale.  This  then. does 
Include  the  “unknowns”.  The  problem  often  arises  where  some 
damage  or  a  manufacturing  defect  beyond  those  that  are  inherent 
occurs  and  an  enalysis  is  required  to  determine  the  suitability 
of  the  chamber  for  use.  This  analysis,  if  limited  to  two- 
dimensions,  may  miss  a  critical  stress  condition  in  the  triaxial 
condition  that  actually  exists. 

Because  of  the  multiplicity  of  possible  failure  modes,  the 
prediction  of  the  ultimate  failure  strength  is  very  difficult 
from  the  theory.  The  original  designs  have  built-in  factors  of 
safety  and  empirical  factors  regarding  the  actual  test  data  of 
similar  designs.  These  factors  cannot  be  directly  applied  to  the 
theoretical  optimization  calculations.  Fundamental  data  will 
have  to  be  obtained  in  order  to  reach  an  optimum  design  condi¬ 
tion.  Crit.*  ~nl  design  optimizations  can  utilize  successive 
approximations  to  create  the  proper  angles  of  windings,  dome 
contours  and  fiber  layers,  as  well  as  resin  content.  The 
routines  to  accomplish  this  can  be  create.'  with  significant 
pay-back  on  the  investment  by  providing  designs  fch't  approach  the 
maximum  possible  performance  and  reliablility  witj.  given  mat¬ 
erials. 
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I  reconend  that  current  deaignera  adapt  the  testing  and 
design  Methods  that  utilize  3-DFE  aa  their  foundation.  The 
hardware  and  software  required  to  accomplish  this  exists.  The 
co^ts  are  justifiable.  Software  to  do  limited  3-DFE  is  currently 
available  for  PC's  that  can  be  on  the  designer's  desktop.  Kore 
advanced  CAD  routines  incorporate  3-dimensional  graphics  that  can 
display  the  grids  and  visually  display  the  stresses  that  a  design 
is  experiencing  to  the  designer/engineer .  These  tools  for  design 
exist  and  should  be  incorporated  into  our  current  design/evalua¬ 
tion  techniques. 


APPENDIX  1 


General  Notea  on  Laminate  Theory  from  "Composite  Materials 
HandbocK" -hel  M.  Schwartz,  McGraw-Hill 

Laainate  analysis  is  a  mathematical  means  of  determining 
elastic  properties  of  the  complete  laainate  from  the  angles  of 
fiber  orientation  and  elastic  moduli  of  the  individual  layers. 
Conversely,  the  method  is  also  used  to  determine  stress  and 
strains  in  each  layer  produced  by  stretching  and  bending  forces 
on  the  complete  laminate. 

The  elastic  properties  of  a  laainate  (Ex,Ey ,Gxy » Yxy>  are 
determined  in  three  steps: 

1)  The  elastic  properties  of  each  layer  are  determined  from 
the  moduli  En,  El2*  and  G12  and  Poisson  ratio  v  12  of  the  layer 
material.  These  properties  are  determined  from  the  Hooke's-law 
relationship  between  stress  and  strain  in  the  principal  material 
directions  (parallel  and  perpendicular  to  the  fiber  direction) . 

2)  The  elastic  properties  of  each  layer  in  the  principal 
laminate  directions  (parallel  and  perpendicular  to  the  axis  of 
the  laminate  structure)  are  determined  from  the  material - 
direction  properties  and  the  angle  between  the  laainate  and  the 
material  directions. 

3)  A  summation  of  the  lamlnate-dlredtion  elastic  properties 
of  all  layers,  taking  into  account  their  relative  positions  with 
respect  to  the  midplare  of  the  complete  laainate,  yields  the 
elastic  properties  of  the  complete  laminate  (Fig.  2) . 

Laminate  analysis  is  also  used  to  determine  stresses  and 
strains  in  each  layer  due  to  forces  on  the  laminate.  Layer 
strains  are  determined  form  laminate  strains  and  rotations, 
taking  into  account  the  layer  position  relative  to  the  middle 
surface  of  the  complete  laminate.  Stresses  in  each  layer  are' 
then  determined  from  layer  strains  and  the  Hooke's-law  relation¬ 
ship.  Thus,  the  stresses  and  strains  at  a  point  in  the.  structure 
can  be  calculated  in  every  direction,  that  is »ax,°y#Txy  and  a  11, 
a22» T12  In  every  ply. 
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V  V 
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J. 


i  on  Linear  Elgstic  Fracture  Mechanics  <LEFM)  Theory 


fro*  “Fracture  Toughness  Testing  and  Its  Applications-Symposiu* 
on  Fracture  Toughness  Testing  and  Its  Applications*  ASTN  Special 
Technical  Publication  No.  381*  American  Society  For  Testing 
and  Materials ,  1964 


The  introduction  of  stresses  in  a  body*  caused  by  a  crack* 
is  a  feature  of  the  strength  of  the  material.  The  stresses  near 
the  crack  tip  may  result  in  the  growth*  or  perpetration  of  the 
cxack  or  flaw.  The  nonlinear  and  plastic  effects  are  neglected 
in  the  gross  features  discussed  in  this  presentation*  but  may  be 
significant  in  a  given  scenario.  The  stress-field  expansions  for 
the  crack  tips  are  applicable  for  the  general  case  of  an  isotro¬ 
pic  elastic  body.  The  three  modes,  opening,  edge-sliding  and 
tearing*  (I*II  and  III  respectively)  are  shown  in  the  figure  (3). 


The  Basic  Modes  of  Crack  Surface  Displacements'. 

Figure  <3>« 


The  surfaces  of  a  crack  dominate  the  distribution  of 
stresses  in  the  vicinity  of  the  crack.  Other  boundaries  are 
remote  and  affect  only  the  intensity  of  the  iweml  stress  field* 
as  does  the  loading  forces  applied  to  the  body.  The  opening 
mode*  I*  is  associated  with  the  local  displacement  where  the 
crack  surfaces  move  directly  apart.  The  edge  sliding  mode*  II* 
is  where  the  displacements  move  over  one  another*  perpendicular 
to  the  leading  edge  of  the  crack.  The  tearing  mode*  III*  is 
where  the  crack  surfaces  slide  parallel  to  the  leading  edge  of 
the  crack. 

Modes  1  and  II  can  be  treated  a*  piane'-extenalonal  problems 
in  the  theory  of  elasticity*  symmetric  and  skew-symmetric, 
respectively.  Mode  III  is  a  pure  torsion  or  shear  problem. 
Using  the  notation  in  figure  <4),  the  resulting  stress  and 
displacement  fields  are  derived: 


Mode  I . 


ox  a  <Ki/(2irr)l/2>CoaO/2Cl-ain<0/2)ain(3©/2)3 
ay  =  <Ki/<2irr)l/2)COa0/2Cl^ain<0/2)ain<30/2)] 
t  xy  =  (Ki/<2iTr)l/2)ain(Q/2)coa<0/2)co3<3  0/2) 

02  =  v(0X  ♦oy > ,  TXy  sTyz  =  0 

u  =  CKi/G)  Cr/<2ir>]  1/2  coa  <0 /2)  * 

Cl -2 v  ♦sin2(0/2>  3 

=  (Ki/G)  Cr/ C2tt)  ]  1/2  ain  (0/2)  * 

C2-2  v  -coa2 (0/2) 3 

w  =  0. 

Mode  II: 

°x  a  -  CKn/<2Tir)l/23aln<  Q/2)  £2*coaC  0/2)coa(3  <^2)3 
oy  *  CKii/(2TTr)l/2]a.'n<0/2)coa<  0/2)coa<3  0/2) 

T  xy  ■  CKu/<2Trr)l/2]  coa  <  0/2)  Cl -ain  C  O'  2)coa(30/2)  3 

°z  B  v  <  °x  *.°y>*  T  xz  *  T  yz  ■  o 
u  «  CKii/G)  Cr/<2ir)  Jl/2  ain<0/2)  * 

C2-2  v ♦co«2<0/2>3 

v  «  <Kn/G>  Cr/<2if)Jl/2  coa<0/2>  • 

1-1*2  v  ♦ain2<9/2)3 

w  »  0. 

Mode  III: 

Txz  *  ~<Km/<2irr) l/2)ain<0 /2) 
ty2  »  <Kni/(2wr)l/2)coa<0/2) 
ox  «  oy  ■  02  •  Txy  •  0 
w  -  <Km/G>  C<2r>/  »ll/2  atn  <  0/2) 
u  •  v  »  0. 
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Coordinates  of  the  Stress  Components  In  the  Crack  Tip. 

Figure  (4) . 

The  parameters,  Kj,  Ku,  and  Km  in  the  equations  are 
stress-intensity  factors  for  the  corresponding  three  types  of 
stress  and  displacement  fields.  These  factors  must  contain  the 
loading  force  magnitudes  linearly  (for  linear  elastic  bodies)  and 
depend  on  the  configuration  of  the  body  including  the  crack 
size.  They,  therefore,  reflect  the  redistribution  of  stress  in  a 
body  due  to  the  introduction  of  a  crack  and  the  force  transmis¬ 
sion  through  the  crack  tip. 

It  has  been  shown  that,  for  the  general  homogeneous  aniso¬ 
tropic  case,  the  crack- tip  stress  fields  and  their  intensity 
factors,  the  complete  analogy  with  the  isotropic  case  is  preserv¬ 
ed.  By  judicious  definition  of  the  anisotropic  stress-intensity 
factors,  they  are  identical  to  those  for  the  isotropic  case. 
However,  for  the  non -homogeneous  anisotropy,  the  resultant  stress 
fields  and  intensity  factors,  as  well  as  stress  distributions, 
may  well  be  different  and  result  in  different  singularities  than 
those  observed  in  the  homogeneous  case. 

The  introduction  of  viscoelasticity  into  the  problem  results 
in  similar  intensity  factors  and  stress  distributions  in  the 
vicinity  of  cracks;  however,  this  causes  these  factors  to  now 
become  variables  as  a  function  of  time. 

The  complete  discussion  and  development,  as  well  as  exten¬ 
sions  to  other  and  general  cases  is  found  in  the  above  credited 
text  and  thi*  presentation  is  merely  s  condensation  of  the 
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principle  governing  equations.  It  is  given  for  completeness  and 
the  readers  reference. 
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